Strain modulation of transport criticality in heterogeneous solids 
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A vast class of disordered conducting-insulating compounds close to the percolation threshold is 
characterized by nonuniversal values of transport critical exponents. The lack of universality im- 
plies that critical indexes may depend on material properties such as the particular microstructure 
or the nature of the constituents, and that in principle they can be influenced by suitable applied 
perturbations leading to important informations about the origin of nonuniversality. Here we show 
that in Ru02-glass composites the nonuniversal exponent can be modulated by an applied mechan- 
ical strain, signaled by a logarithmic divergence of the piezoresistive response at the percolation 
threshold. We interpret this phenomenon as being due to a tunneling-distance dependence of the 
transport exponent, supporting therefore a theory of transport normniversality proposed some years 
ago. 
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PACS numbers: 72.20.Fr, 64.60.Fr, 72.60.+g 



Transport properties of disordered insulator-conductor 
composites are characterized by the existence of a perco- 
lation critical value Xc of the conducting phase volume 
concentration x below which the system becomes an in- 
sulator. As a; — sic —> the resistivity p of the composite 
diverges by following a power-law behavior of the form: 



p ~ pq{x - XcY 



(1) 



where po is a material-dependent prefactor and t is the 
transport critical exponent P, According to the stan- 
dard theory of transport in isotropic percolating materi- 
als, po and Xc depend on microscopic details such as the 
microstructure and the mean inter-grain junction resis- 
tivity while, on the contrary, the exponent t is material 
independent and takes on the value t ~ 2.0 for three di- 
mensional systems 0, This universal value is indeed 
observed in a great variety of composites 0, 0, 0, |^ and 
it is confirmed to a rather high accuracy by calculations 
on random resistor networks models 0- 

Universality of DC transport is understood as being 
due to the unimportance of microscopic details for the 
flow of current compared to the topological properties of 
the percolating backbone fl)!. This concept applies also to 
other transport-percolation phenomena such as for exam- 
ple proton conductivity in biological materials and lami- 
nar fluid flows through porous media 0]. Yet, a vast class 
of disordered composites deviates from universality by 
displaying values of t larger than 2.00, |^ and as high 
as t ~ 7.0^3 or even more(ll!|. The lack of universal- 
ity indicates that the exponent t acquires an additional 
dependence upon some microscopic properties and vari- 
ables, whose identification would permit to unveil the 
mechanisms leading to t ^ 2.0. A practical route to 
achieve this goal is to apply an external perturbation to 
a nonuniversal composite with the hope of affecting the 
same microscopic variables which govern nonuniversality 
In this case, the applied perturbation would change the 
value of the critical exponent t. 



We have applied this idea to Ru02-glass composites 
prepared as thick-film resistors (TFRs) in our laboratory. 
Ru02 as well as other ruthenate (Pb2Ru206, Bi2Ru207) 
or oxide (Ir02) TFRs are ideal systems for our study 
because they display universal {t ~ 2.0) or nonuniversal 
(t > 2.0 6] up to about t ~ T.oflfij) behavior of trans- 
port depending on the microstructure and fabrication 
procedures, with a good reproducibility and high sta- 
bility. In these composites, the conducting particles are 
separated by nanometer-thick film of glass[l3 and the 
dominant mechanism of transgort is via quantum tun- 
nehng through the glass films Here we show 

that when transport is nonuniversal, the critical expo- 
nents of our samples can be modified by an applied me- 
chanical strain e, providing evidence that the microscopic 
variable governing nonuniversality is the mean tunneling 
distance a between adjacent conducting grains. 

Our samples were prepared by mixing two series 
of Ru02 powders with 40 nm and 400 nm grain 
sizes, respectively, with a lead-borosilicate glass powder 
[PbO(75% wt)-B2O3(10% wt)-Si02 (15% wt)] of 1 - 5 
/im grain sizes together with a vehicle of terpineol and 
ethyl cellulose. 2% of AI2O3 were added to the pastes 
in order to avoid crystallization. The pastes were screen 
printed on AI2O3 substrates with gold electrical contacts 
and fired for 15 minutes at temperatures Tf (see Tabled 
higher then the temperature of fusion of the glass (about 
500°C). The resulting films were about 10 pm thick and 
several resistivity measurements were taken over eight 
different samples for each Ru02 volume fraction value. 

In Fig^we report the room temperature resistivity p 
measured for four different series of TFRs (see Tablelill as 
functions of the RUO2 volume concentration a;. As shown 
in Fig^a), resistivity diverges at rather small values of a;, 
as expected when the mean grain size of the conducting 
phase (40 nm and 400 nm) is smaller than that of the 
glass (1-5 pm)^^. The same data are re-plotted in the 
In-ln plot of FigQJb) together with the corresponding fits 
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FIG. 1: (a): resistivity p as a function of Ru02 volume con- 
centration X for four different series of TFRs. The solid lines 
are a guide to the eyes, (b) In-ln plot of the same data of (a) 
with fits to Eq.10 shown by solid lines. The dashed line has 
slope t = 2 corresponding to universal behavior of transport. 
The prefactor po, critical concentration Xc and transport ex- 
ponent t values obtained by the fits are reported in Table 

HI 



to Eq.(^ (solid lines) and the best-fit parameters poi 
and t are reported in Table ^. As it is clearly shown, our 
resistivity data follow the power law behavior of Eq.QJ 
with exponent t close to the universal value t ~ 2.0 for 
the Al series (t = 2.15 ± 0.06) or markedly nonuniversal 
as for the A2 case which displays t = 3.84 ±0.06. The Bl 
and B2 series have virtually equal values of i (i ~ 3.16) 
falling in between those of the Al and A2 series. 

The effect of an applied strain e on transport is ob- 
tained by recording the piezoresistive response, i. e. the 
relative resistivity change ^p/ p upon applied mechani- 
cal strain deformation, by cantilever bar measurements. 
The Ru02-glass pastes were screen printed on AI2O3 can- 
tilever bars 51 mm long, 6=5 mm large, and h — 0.63 
mm thick. The cantilever was clamped at one end and 
different weights were applied at the opposite end. The 
resulting substrate strain e along the main cantilever axis 
can be deduced from the relation e — 6Mgd/ (Ebh^), 
where d = 27.8 mm is the distance between the resistor 
and the point of applied force, E = 332.6 GPa is the 
reduced AI2O3 Young modulus, g is the gravitational ac- 
celeration and M is the value of the applied weight. As 
shown in RefsEl, 113, 

from measurements of the longi- 
tudinal (Ap||/p) and transverse {Ap±/p) piezoresistive 
responses obtained by measuring the voltage drops along 
and perpendicular to the main cantilever axis, respec- 
tively, it is possible to extract the isotropic resistive varia- 
tion defined as Ap/p = {Ap^^/p+2Ap±/p)/3. For the val- 
ues of Ru02 content used in this work, Ap||/p ~ A p 1 /p 
confirming the vicinity to a percolation threshold j 1 611 1 7| . 



FIG. 2: (a) relative variation of resistivity as a function of 
applied strain e in cantilever bar measurements of the A2 
series for different contents x of RUO2. The solid lines are 
linear fits to the data, (b) isotropic piezoresistivity response 
r — d\n{p)/de obtained from the slopes of (a) as a function 
of the volume concentration x for the TFRs series of Fig0 



In FigI2ta) we report Ap/p for various Ru02 con- 
centrations X of the A2 series as a function of applied 
strain e. In the whole range of applied strains, Ap/ p 
changes linearly with e permitting to extract rather accu- 
rate piezoresistive factors T — dln{p)/de from the slopes 
of the linear fits of Ap/p vs e. The so obtained F val- 
ues as a function of x are plotted in Fig|2Ib) for the four 
TFRs series of Fig^ With the exception of the Al series, 
F displays a strong dependence upon x and tends to di- 
verge as X approaches to the same critical concentrations 
Xc at which p goes to infinity. 

The divergence of F at Xc has been reported some years 
agoJSjl, but here we show that it is caused by a strain 
modulation of the critical exponent t. To see this, let us 
consider Eq.Q and differentiate it with respect to e: 



de de 



tin 



(2) 



Since Fq = d\n{pQ) / de is a constant independent of a;, all 
the X dependence of F must come from the last term of 
Eq.(j21l. However, it must be noted that a; is a measure 
of the concentration p of intergrain junctions with finite 
resistances present in the sample 0, 13] ■ In our measure- 
ments, the applied strain is only of order 10~^, which 
is so small that it can eventually change the value of 
the tunneling junction resistances but cannot modify the 
concentration p of junctions. So dx/de must be zero and 
by the same token also dxc/de — 0. Hence the only way 
to have a x dependence of F is to allow the transport ex- 
ponent t to have non vanishing derivative, which reduces 



3 



TABLE I: Label legend of the various sample used in this work with fitting parameters of Eas. lll3ll 



Label 


Ru02 grain size 


firing temperature T/ 




ln(po/nm) 


t 


To 


dt/de 


Al 


400nm 


525° C 


0.0745 


-11.1 ±0.3 


2.15 ±0.06 


5.5 ± 1.5 


-0.2 ±0.4 


A2 


400nm 


600° C 


0.0670 


-14.2 ±0.2 


3.84 ± 0.06 


-8.8 ± 1.6 


5.4 ±0.5 


Bl 


40nm 


550° C 


0.0626 


-14.3 ±0.5 


3.17 ±0.16 


-15.3 ±3 


8.7 ±0.9 


B2 


40nm 


600° C 


0.0525 


-13.7 ±0.7 


3.15 ±0.17 


-19.3 ± 2.4 


11.0 ±0.7 



Eq.© to: 



dt f 1 
r = To + — In 

ae \x — Xr 



(3) 



In Figl^l wc plot the piezoresistive data of Fig|2Ib) as 
a function of ln(a; — Xc) with the same values of the criti- 
cal concentrations Xc extracted from the resistivity data. 
Our data fit well with Eq.© (solid lines and last two 
columns of Table I) and indicate that T indeed diverges 
logarithmically 3tS X — * Xq for the series A2, Bl, and B2. 
An analysis of previously published data'lS] confirms this 
conclusion 19]. It is also remarkable that the series Al 
has vanishing dt/de within experimental errors, in per- 
fect agreement with the universal value of its exponent 
{t ~ 2.15 ± 0.06). Wc conclude therefore that when the 
transport exponents of our Ru02 TFRs are nonuniversal, 
they can be modulated by an applied strain. 
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FIG. 3: Piezoresistive factor as a function of ln(a; — Xc) and 
fits (solid lines) to Eq.Q. The fit parameters dt/de and To 
are reported in Table I. Only sample Al (solid squares) has 
quasi-zero slope indicating that dt/de = 0. 



Let us discuss now what this finding implies in terms 
of the microscopic origin of nonuniversality. The main 
effect of the applied strain on the microscopic proper- 
ties of Ru02-glass composites is that of changing the 
mean intergrain tunneling distance a, a a(l ± e), 
leading to a variation of the microscopic tunneling re- 
sistances. Hence, in order to have dt/de ^ 0, the trans- 
port exponent itself must depend on a. A scenario of 
this kind has been proposed by Balberg a few years ago 
in his tunneling-percolation theory of nonuniversality in 
carbon-black-polymer composites 20] . According to this 
theory, when the distribution function of the tunneling 
distance d between two neighboring grains decays with d 
much slower than the tunneling decay exp(— 2d/^), where 
^ is the localization length, then the distribution function 
h{r) of the intergrain tunneling resistances r develops a 
power-law tail such that h[r) cx r"~^ for r — > oo where 
a — 1 — ^/(2a). It is well known |2l| that distribution 
functions of this kind lead to an explicit dependence of t 
upon the parameter a 22, 23] : 



t 



to 

V - 



1 

l-a 



if + < ^0 

if V+-r^>tQ 



(4) 



where t^ ~ 2.0 is the universal transport exponent and 
~ 0.88 is the 3D correlation-length exponent. In view 
of Eq.|01), the term which multiplies the logarithm in 
Eq.Q is dt/de = 1/(1 - a) 2a/^ > when t is given 
by the second line of Eq.Q) or dt/de = when t — to. 
However, we note that in Ru02-based TFRs the large dif- 
ference between the bulk modulus B values of Ru02 and 
of the glass (BruQs - 270GPa and Bgiass =i 40 - 80GPa) 
leads to local strain variations resulting basically in a 
strong e amplification within the softer phase (glass) 
through which tunneling occurs. Hence, the dt/de val- 
ues of A2, Bl and B2 reported in Table I are not simply 
equal to 2a/^ but incorporate also the amplification effect 
of the strain heterogeneity. 

An interesting feature of our piezoresistive data is the 
change of sign of Fq which from positive for the uni- 
versal serie Al becomes negative for the nonuniversal 
ones (A2,B1 and B2, see TablePl. Although being quite 
surprising at first sight, this result can be natu rally ex- 
plained within the tunneling-percolation model|20|. In 
fact a tensile strain (e > 0) enhances the intergrain tun- 
neling resistances leading to an overall enhancement of 
the sample resistivity, so that F = d\Ti[p)/de must be 
positive. Hence, from Fq.©, Fq > when dt/de = 
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while when dt/de > Fq does not need to be positive to 
ensure F > 0. Quite remarkably, an effective medium cal- 
culation with the tunneling-percolation distribution func- 
tion h{r) oc r"~^ leads to Fq < as long as t depends 
upon a as in the second line of Eq.QjlJj. This is in 
agreement with our experimental results (see Table I). 

Let us now turn to discuss briefly how other theories of 
nonuniversality fit our findings. First, the random void 
model (RVM) of continuum percolation [i^l can hardly 
be apphed to the Al and A2 series, since the Ru02 
particles are not microscopically smaller than the glass 
grains, so that the conducting phase cannot be regarded 
as a continuum. In addition, the RVM predicts that 
t follows Eq.|01) with a = Orvm being a purely topo- 
logical quantity defined by the distribution of the nar- 
row necks bounded by three interpenetrating insulating 
spheres. Hence, aRVM is insensitive to an applied strain 
so that dt/de — and the piezoresistivity F cannot de- 
velop a logarithmic divergence. We exclude also the pos- 
sibility that our nonuniversal samples are in the mean- 
field regime for which t = tmi = 3 [2^. In fact the A2 
series has a critical exponent clearly larger than and, 
more importantly, dtrai/de is expected to be strictly zero. 

In summary, we have shown that the piezoresistive re- 
sponse of disordered Ru02-glass composites has a log- 
arithmic divergence at the percolation threshold when 
DC transport is nonuniversal. A coherent interpreta- 
tion of this result calls into play a mean tunneling dis- 
tance dependence of the resistivity exponent, in agree- 
ment with a tunneling-percolation origin of nonuniver- 
sality proposed some time ago 20]. Such mechanism of 
universality breakdown could apply also to other materi- 
als for which transport is governed by tunneling such as 
carbon-black-polymer composites, and experiments on 
their piezoresistive response could confirm such conjec- 
ture. 

Before concluding, it is worth to point out that in addi- 
tion to the importance of our results concerning the origin 
of transport nonuniversality in tunneling-percolation ma- 
terials, the reported logarithmic divergence of F has also 
consequences for applications of Ru-based TFRS in force 
and pressure sensors. In fact, the efficiency of piezoresis- 
tive sensors is measured by the magnitude of F. However 
the logarithmic increase of F as x — > is not sufficient 
to overcome the power-law divergence of the resistance 
fiuctuations. This poses serious limitations regarding 
the precision of pressure/force detection. An alternative 
route would be to make use of the amplification effects 
induced by the local elastic heterogeneities by fabricating 



Ru02-based TFRs with softer insulating phases like low 
melting temperature glasses or even polymeric hosts. 
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